Heavy-light meson spectrum with and without NRQCD 
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Results for the spectrum of S and P-wave charmed mesons are obtained in the quenched approximation from 
a tadpole-improved anisotropic gauge field action and a D234 quark action. This is compared to the spectrum 
obtained from an NRQCD charm quark and a D234 light antiquark. NRQCD results for bottom mesons are also 
discussed. 



1. MOTIVATION 

If the charm quark is sufficiently heavy, then 
lattice NRQCD is an efficient computational 
method for obtaining the spectrum of mesons 
with a single charm quark. Recent studies of 
the S and P-wave meson spectrum have deter- 
mined the contributions from first, second and 
third order in the inverse charm mass expan- 
sion. [|U|] Second and third order contributions 
were found to be statistically insignificant for the 
P-wave mesons, but the S-wave spin splittings 
(D* — D and D* —D s ) acquire substantial correc- 
tions at both second and third order. It is diffi- 
cult to make a firm statement about convergence 
of the expansion from these findings alone. 

In the present work, the existing results are 
compared to a new study, where the NRQCD 
charm quark is replaced by a relativistic charm 
quark using a D234 fermion action. @ The new 
simulations use the same set of gauge field con- 
figurations and the same light antiquark propa- 
gators as were used in Ref. ^| , so the only differ- 
ences are the heavy quark action (NRQCD versus 
D234) and the operators used to create or destroy 
the various meson states. 

The S-wave masses have been discussed many 
times in the literature (see Refs. Jlj^Q and refer- 
ences therein), but the P- waves have only recently 
been extracted from lattice data: the limit of an 
infinitely heavy quark was studied by Michael and 
Peisa||, a subset of the P-wave charmed masses 
were found by Boyle using the clover action 



(he also extrapolated to the bottom region) 
and two groups have used NRQCD to study the 
charmed and bottom meson spectra[^],[|,[7]||. 

2. LATTICE CHOICES AND METHOD 

The gauge fields of Ref. 0, also used for the 
present study, are on 10 3 x 30 anisotropic lattices 
with (a s /at )bare = 2. The gauge action is classi- 
cally improved by including rectangular plaque- 
ttes, and tadpole improved using the mean link in 
Landau gauge. The coupling is fixed at (3 = 2.1, 
and the renormalized anisotropy is found to be 
a s /a t = 1.96(2). 

The light antiquark fields are also taken from 
Ref. p| . They are derived from a D234 ac- 
tion Upwhich removes leading and next-to- leading 
classical errors and is tadpole improved. Two 
hopping parameters (k — 0.23 and 0.24) are 
used, corresponding to "m w /m p " = 0.815(3) and 
0.517(8). The physical p meson mass leads to 
a t — 0.1075(23) fm, and the physical kaon mass 
implies K strangc = 0.2356(3). Dirichlet time 
boundaries are used for the light antiquark. 

In Ref. |H the heavy quark was described by 
NRQCD, and new results are presented here 
where a heavy D234 quark is used with k = 0.182, 
which is in the vicinity of charm. By interpolat- 
ing the NRQCD results to this same heavy quark 
mass, a direct comparison of NRQCD and D234 
results can be made. 

With and without NRQCD, heavy- light me- 
son masses come from 2000 quenched configura- 
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tions, and statistical uncertainties are obtained 
from 5000 bootstrap ensembles. The heavy-light 
meson creation operators are 

J2QH^(x)[l + c s AW(x)] n °q(x) (1) 

X 

for NRQCD, and 

$(x)n(x) [1 + c s F(x)] ns ip{x) (2) 
for D234, where 

3 

F(x) = £ [Ui(*)* x ,y-i + Uj(x - i)S x ^ (3) 
t=i 

and A^ 2 ) is the lattice Laplacian. 

The quantum numbers are dictated by fl as 
presented in Table |j} At the source, the parame- 
ters (c s ,n s ) are set to (0.15,10) for NRQCD and 
(1,15) for D234. The sink operator is local in all 
simulations. 

3. S-WAVE HEAVY-LIGHT MASSES 

The S-wave hyperfine splitting for D234 is 
shown in Fig. |l| in lattice units, for the two avail- 
able light antiquark hopping parameters. Also 
shown are the NRQCD results from Ref. Q after 
interpolation of the 1 So kinetic mass to match the 
D234 value of the 1 S'o meson mass. This interpo- 
lation was done independently at each order in 
the NRQCD expansion, so the 1 S'o mass main- 
tains a common value in each case. From Fig. |l|, 
lowest order NRQCD is seen to give smaller hy- 
perfine splittings than D234. The addition of 
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Figure 1. Comparison of D234 and NRQCD 
for the S-wave spin splitting, in lattice units. 
Solid(open) symbols use k = 0.23(0.24) for the 
light antiquark. 

0(1/M 2 ) terms tends to increase the NRQCD 
results, but 0(1/M 3 ) contributions are larger in 
magnitude than the preceding order so the status 
of the NRQCD expansion for this observable is 
unclear. 

The NRQCD results are shown in physical 
units in Fig. ||. Whereas Fig. [j] compared the 
different NRQCD orders by choosing a common 
physical 1 S'o mass, Fig. || shows the NRQCD ex- 
pansion when the bare charm quark mass is held 
fixed. Neither plot gives a compelling defense 
for a convergent 1 /M c h arm expansion, but conver- 
gence is not disproven by such plots either. The 
D s — D splitting receives insignificant corrections 
from second and third orders, but corrections to 
the spin splittings are significant. Experimentally 
the D s — D + splitting is 104 MeV, and therefore 
agrees with the NRQCD determination, but the 
spin splittings are D* + - D + = 141 MeV and 
D* s -D s = 144 MeV. With or without NRQCD, 
the quenched prediction is smaller than experi- 
ment. 

Ref. |Q also provides results for the bottom 
mesons, and finds that the S-wave bottom masses 
are completely dominated by leading order in 
NRQCD so convergence is not disputed. As for 
charm, the spin splittings are smaller than exper- 
iment. 
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Figure 2. NRQCD S-wave mass splittings in 
physical units. 



4. P-WAVE HEAVY-LIGHT MASSES 

In Ref. nonleading corrections to P-wave 
heavy-light masses were found to be insignifi- 
cant for both bottom and charm mesons. This 
leads one to expect agreement between D234 and 
NRQCD determinations of these P-wave masses. 
Fig. U shows all four of the P-wave masses in lat- 
tice units for both of the available light antiquark 
hopping parameters. D234 results are new, and 
the NRQCD results are interpolations of the data 
from Ref. Q so that the 1 Sq kinetic mass matches 
the D234 1 Sq mass. Only statistical uncertainties 
are shown. 

According to Fig. ||, the P-wave masses do not 
depend upon whether the charm quark is de- 
scribed by NRQCD or D234. One rather dra- 
matic exception seems occur for the 1 P\ with 
the heavier antiquark, but this systematic error 
is understood: the D234 calculation has no true 
plateau in this case. The plateau-finding algo- 
rithm of Ref. |^| has been carried over to the 
present work, and it chooses the best plateau by 
maximizing the quality factor, 

^_ r(7V/2-l, X 2 /2) 



r(7v/2 - i,o) 



where 
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Figure 3. Comparison of D234 and NRQCD for 
the S-P splittings, in lattice units. Solid(open) 
symbols use k — 0.23(0.24) for the light anti- 
quark. 



and N is the number of timesteps in the proposed 
plateau region. The D234 1 P\ meson becomes 
quite noisy at early Euclidean times such that no 
clear plateau is evident, so the method of max- 
imum Q chooses a "plateau" which begins too 
near the source. This leads to the erroneously 
large mass shown in Fig. |j. Except for this one 
correlator, Fig. |^ presents a clear agreement be- 
tween results with and without NRQCD. 

The 3 Po and 3 P2 masses for charm and bottom 
are displayed in Figs. ^ and || from the NRQCD 
results of Ref. f§], along with NRQCD results 
from Ref. J§| , the clover results of Ref. || , and a 
variety of model calculations || . A disagreement 
between the two lattice NRQCD calculations for 
£?2 — Bq is evident; Ref. [g| notes that no choice 
of plateau region allows those data to attain the 
large splitting of Ref. The lattice results of 
Ref. Q cannot resolve the discrepancy. The non- 
lattice models are more easily distinguished from 
one another for charm than for bottom, and the 
lattice result clearly favours a small D% — Dq split- 
ting (substantially less than 100 MeV). 
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Figure 4. Splittings between the 3 P2 and 3 Pq 
masses for the D** and B** systems. Open(solid) 
symbols involve an s(u,d) quark. 



Figure 5. S-P splittings for heavy-light mesons. 
Open(solid) symbols involve an s(u,d) quark. 
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